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ABSTRACT 

Despite over 30 years of study, the mass-area relationship within and among clouds is still poorly 
understood both observationally and theoretically. Modern extinction datasets should have sufficient 
resolution and dynamic range to characterize this relationship for nearby molecular clouds, although 
recent papers using extinction data seem to yield different interpretations regarding the nature and 
universality of this aspect of cloud structure. In this paper we try to unify these various results 
and interpretations by accounting for the different ways cloud properties are measured and analyzed. 
We interpret the mass-area relationship in terms of the column density distribution function and its 
possible variation within and among clouds. We quantitatively characterize regional variations in the 
column density PDF. We show that structures both within and among clouds possess the same degree 
of "universality" , in that their PDF means do not systematically scale with structure size. Because 
of this, mass scales linearly with area. 

Subject headings: ISM: structure, ISM: clouds, stars:formation 



1. INTRODUCTION 

The hierarchical spatial-kinematic structure of molecu- 
lar clouds provides clues about the processes that dictate 
cloud evolution, and sets the initial conditions for star 
formation. Cloud structure is influenced by several pro- 
cesses - magneto-hydrodynamic turbulence, gravity, stel- 
lar feedback, internal and external radiation fields, etc. 
- and is extremely complex, exhibiting features across a 
wide range of scales. Absent a framework for describing 
cloud structure in its entirety, researchers have applied a 
bevy of empirical techniques to characterize observed and 
simulated data; examples include the distribution func- 
tions for density, column density, velocity and mass, size- 
linewidth relationship, structure function, delta- variance 
relationship, principal co mponents analysis, and spe ctral 
correlation function (see lElmegreen fc Scald 12001 and 
references therein). The relationships and interdepen- 
dences that exist between these various descriptors are 
rarely obvious. 

Larson's scaling relationships ()Larsonlll981|) are among 
the most versatile empirical characterizations of molec- 
ular cloud structure. The third of these relationships 
(which we will refer to as L3) states that the average vol- 
ume density of a molecular cloud scales inversely with its 
size or, alternately, that a cloud's mass is proportional to 
its area. In the thirty years since its publication, Larson's 
mass-area relationship has been subject to numerous re- 
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investigations. Table [TJ lists several published measure- 
ments of mass-area relationships, assuming a functional 
form (M/Mq) — 7 (A/pc 2 ) 13 . These are all measure- 
ments of Galactic cores, clouds, a nd cloud complexes 
(for extragalactic m easurements, see lBolatto et alj|2008t 
IHughes et alj 12010). There is substantial scatter among 
these values, especially in the constant of proportionality 
7- 

Much of the uncertainty about L3 (and the scatter 
in Table [J) is due to the heterogeneous ways in which 
cloud properties are measured. The conversion from in- 
tensity to mass is non-trivial for most tracers; dust emis- 
sion is influenced by uncertain temperature and emissiv- 
ity variations, while molecular line emission is affected 
by op acity, chemistry, and changing excitation condi- 
tions (|Goodman et alj |2009). Several authors have sug- 
gested that L3 is an observational artifact, owing to the 
fact that individual surveys are sensitive to a relatively 
narrow range of surface brightnesses, artificially driving 
measurements towards the M oc A relation ship (jKegell 
119891 : iBallesteros-Paredes fc Mac Low! 120021 ). Further- 
more, clouds do not have sharp edges that define a natu- 
ral size; researchers usually use specific contours or mo- 
ment measurements to describe a cloud's area. Finally, 
it is difficult to reliably segregate clouds located in the 
Galactic mid-plane. All of these factors will affect the 
mass-area relationship. 

Larson's relationships are al so studied in d isparate re- 
gions. Indeed, the sample in iLarsonl (|1981f ) comprises 
small substructures within clouds, individual clouds, and 
large cloud complexes. It is important to distinguish 
cloud-to-cloud comparisons from comparisons of struc- 
tures within clouds, as these two regimes are governed 
by different processes and may yield different scaling re- 
lationships. 

Lombardi, Alves and Lada ((2JDlfj, hereafter LAL10) 
addressed many of these issues by measuring cloud sizes 
and masses in a syste matic way. Using 2MASS-derived 
(jSkrutskie et al.1 120061 ) extinction maps (which do not 
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suffer the same dynamic range limitations or calibra- 
tion uncertainties inherent to dust and gas emission 
data; iGoodman et al.l [2009) , they measured the mass 
and area of 11 nearby clouds. All of these clouds are 
far enough out of the mid-plane that background con- 
fusion is low. They experimented with different cloud 
boundaries, using the extinction contours at Ak = 
(0.1,0.2,0.5,1.0,1.5). Comparing different clouds at a 
constant extinction threshold, they found that the data 
obey M oc A to excellent approximation, although the 
constant of proportionality changes with each threshold 
(Figure QJ,) . Their results provide strong evidence that 
the correlation between the areas and masses of clouds is 
a real phenomenon, and not an observational bias. They 
also demonstrate that a lognormal column density dis- 
tribution (PDF) can reproduce this relationship, as long 
as the cloud-to-cloud variation in the lognormal width 
parameter (er) is not too large (they report a scatter of 
~ 20%). From their results, they conclude that molecu- 
lar cloud structure is "universal", in the sense that dif- 
ferent clouds possess similar column density PDFs, and 
as a result of this obey M oc A. 

Kauffmann ct al. (2010a. b) have also recently studied 
the mass-area relationship of nearby clouds using extinc- 
tion data. Like LAL10, these studies defined cloud sub- 
structures by examining regions with extinctions above 
a given threshold. However, Kauffmann et al. ana- 
lyzed each closed contour separately; Lombardi, Alves 
and Lada merged all regions above a given contour level 
when measuring cloud masses and sizes (see Figure [5]). 
Superficially, the mass-size relationships that Kauffman 
et al. derived within clouds seem to exhibit more het- 
erogeneity than the relationships that LAL10 derived be- 
tween clouds (Figure [Tp). This raises a question: to 
what extent is the "universality" in cloud structure em- 
phasized by LAL10 (and attributed to the global simi- 
larity of cloud column density PDFs) challenged by the 
internal structural heterogeneity suggested by the work 
by Kauffmann et al.? 

Answering this question requires more carefully com- 
paring the methodological differences between these 
studies, and is the focus of this paper. We emphasize the 
distinction between the mass-size relationship between 
and within clouds, and interpret both in the context of 
the column density PDF and its variations. We first 
consider the inter-cloud relationship, and provide a few 
important critiques to LAL10 regarding the conditions 
under which a collection of column density PDFs will 
yield L3. These critiques lead to a more specific defini- 
tion of "universal" cloud structure. We then consider to 
the intra-cloud mass-area relationship. This relationship 
can also be connected to the column density PDF, and in 
particular can diagnose variations in the PDF within the 
cloud. We demonstrate that the substructures within 
Perseus possess the same quality of universality as do 
structures across clouds. 

2. DATA AND ANALYSIS 

In this paper we use the same extinction data pre- 
sented by LAL10. These maps are produced using 
the NICEST technique of measuring the reddening of 
background star s due to foreground molecular clouds 
(iLombardil 120091) The maps are based on 2MASS data 
Skrutskie et al.1 120061) . have an effective angular reso- 



lution of 1 — 3', and a noise of Ak ~ 0.06 (N ~ 
5 x 10 20 cm -2 ) - comparable to the median column den- 
sity within mol e cular clouds. We refer the reader to 
lLombardi et al.1 (|2010l) for more details about the data 
reduction process. 

We define cloud structures via contours in the extinc- 
tion maps. The area of each closed contour is obtained 
directly by multiplying the surface area of each pixel by 
the number of pixels within the region. Likewise, the 
mass is obtained by integrating the extinction within 
the region, and using the re lationship 1 mag Ak = 
180M Q pc- 2 (jRieke et al.lll985D . Many authors further 
define a cloud's radius via R = \J A/ir. However, to avoid 
confusion when dealing with significantly non-circular re- 
gions, we restrict our focus to area measurements only. 

As mentioned above, LAL10 merge together all pixels 
above a given threshold when measuring mass and size 
- that is, at a given threshold, each cloud is described 
by a si ngle mass and size (Figure [5]). Kauffman n et all 
(2010aj), on the other hand, treat each closed contour 
separately, and thus are sensitive to potential variations 
within a cloud. 

Several mass-area relationships can be measured given 
a collection of clouds and their substructures. These are 
depicted schematically in Figure [31 where each point cor- 
responds to a mass and size measurement of a single 
closed contour. The shape of the point depicts which 
cloud each measurement comes from, and the color de- 
picts the contour value used. Each relationship is defined 
by whether or not the comparison spans several clouds, 
and whether structures are defined at fixed or variable 
thresholds. Table [T] lists the comparison types used in 
the literature. We restrict our focus to comparisons of 
structures defined at a common extincti on threshold (i.e. 
the T ypes 2 and 4 in the Figure). As lLombardi et all 
(2010) demonstrated, the choice of extinction thresh- 
old affects the offset in the mass-area relationship, since 
the same region defined using a lower extinction thresh- 
old will contain more a r ea an d mas s. Furthermore, 
as bot h lLombardi et al.l ([20101 ) and iKauffman n et al.1 
(|2010al ) note, comparing structures defined using differ- 
ent thresholds will bias the relationship, since structures 
defined at low thresholds sample systematically lower 
surface-density material. 

3. THE INTER-CLOUD MASS-AREA 
RELATIONSHIP 

When cloud structures are identified and measured us- 
ing extinction contours as described above, their masses 
and areas are related by 



M = A (N) A 



(1) 



where (N) is the mean column density inside that con- 
tour, and A is a constant of proportionality that converts 
column density (particles per area) to surface density 
(mass per area). We stress that this equation does not 
necessarily imply that M ex A, since the quantity (N) 
can vary from region to region. As LAL10 note, the 
mean column density can also be expressed as an inte- 
gral of the column density distribution function P(x): 



(N) 



xP(x)dx 



(2) 



N 
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Here, Nq is the column density threshold used to define 
the cloud boundary, and P(x) is the normalized column 
density distribution function - the units of the function 
are inverse column density, and its integral over all col- 
umn densities is one. Importantly, P(x) describes only 
the distribution of pixels within the region of interest - 
this dis tribution can, in princ iple, vary from region to 
region. lLombardi et al.l (|2010f ) circumvent any possible 
regional PDF variation by treating all pixels above Nq 
as a single entity, even if this threshold divides the cloud 
into several disconnected regions (Figure [2]). In this way, 
the distribution function P(x) used in Equation [2] is al- 
ways the same for a single cloud. 

It is clear from Equation [T] that Larson's mass-area re- 
lationship (i.e., M oc A) is approximately satisfied when- 
ever (N) is approximately constant for the objects stud- 
ied. Alternately, if a common threshold No is used to 
define all objects in question, L3 is approximately sat- 
isfied whenever P(x) is approximately constant across 
regions. 

The clouds presented in lLombardi et al.l (|2010|) exhibit 
excellent agreement with L3. The authors further show 
that, if clouds possess log-normal column density PDFs, 
then the parameters describing the location and width of 
this distribution must be similar across regions to repro- 
duce this result. The log- normal fits to these PDFs ex- 
hibit 60% and 20% scatter in the location (/i) and width 
(a) parameters, respectively. It is important to empha- 
size that the lognormal distribution is not the only distri- 
bution capable of reproducing L3 (indeed, any P(x) can). 
Generalizing the analysis in LAL10 to arbitrary PDFs 
reveals interesting connections between the column den- 
sity distribution and mass-area relationship, and more 
precisely defines under what conditions L3 holds. 

Consider a collection of mass and area measurements 
(Mi, Ai) for q regions, using a common threshold Nq, and 
a power-law model: 

log Mi = a + b x log Ai (3) 

We seek the conditions under which the data are well 
described by L3 - i.e., Equation [3] with b — 1. Note that 
we can re-arrange Equation [I] as follows: 

log M, = log A + log (N) + 1 x log Ai + a (4) 

where log (N) is the sample mean of the individual 
log (N) terms, and = log (N) i — \og{N). Equations 
[3] and |4] have the same form, and differ by a residual 
term e^. If the ti terms are uncorrelatecQ, then Equa- 
tion [3] is an appropriate model, and a least-squares fit 
is expected to recover a = log A + log (N) and b = 1 to 
within some statistical uncertainty. This restriction on e 
mandates that log (-/V) must not correlate with A. Equiv- 
alently, the mean of the column density PDF must not 
scale with cloud size for Equation [3] to be an appropriate 
model. If this is the c ase, then the statis tical uncertainty 
on b will be given by (jPress et al.ll2007l section 15.2) 



7 In a formal least-squares analysis, the t{ terms must be inde- 
pendent and drawn from a zero-mean Gaussian distribution. We 
have relaxed the requirement that their distribution be Gaussian, 
and thus treat Equation \5\ as an approximation. 



Here q is the sample size, and c>f og ^ and &\ og A are 
the sample variances of the log (N) i and log Ai measure- 
ments, respectively. Thus, the precision at which L3 is 
recovered will be high when the dispersion of PDF means 
is small compared to the dispersion in areas. In summary, 
the conditions sufficient for a sample of mass and area 
measurements to follow L3 with low scatter are: 

(N) A 

°"lo g (JV> <<Z °f og A ( 6 ) 

We use Equation [5] to determine when L3 is satis- 
fied, and use these criteria as a more precise opera- 
tional definition for cloud structure universality. These 
re sults provide thr e e crit iques to the analysis presented 
in lLombardi et all (|2010t ). First, as noted above, L3 is 
not specific to the log-normality of cloud column density 
distributions. Second, the slope of the mass-area rela- 
tionship is influenced not only by the amount of total 
variation in the PDF, but also by any systematic scaling 
of the PDF with cloud size. Finally, in addition to the 
variation in PDF, both the sample size and dispersion 
in sizes influence the expected scatter about the M oc A 
relationship. 

These critiques are illustrated in Figures |4][7] Each fig- 
ure shows a collection of column density PDFs, the size 
and mean column density of each region, and the mass- 
area relatio nship. For reference, Fi gure |U shows the data 
presented in[Lpmbardi et al.l (j2010f ) . There are mild vari- 
ations in the PDF for each cloud (left), but these varia- 
tions introduce only minor scatter in the measurements 
of (N), which are uncorrelated with size (right). Con- 
sequently, the clouds obey L3 (bottom). Figure [5] shows 
an example of hypothetical structures with pathologi- 
cal, non-lognormal PDFs. These PDFs are also similar 
enough that (N) is roughly constant, and these clouds 
also obey L3. Figure |6] presents a modified version of the 
data in Figure SJ Here, the measurements of A are shuf- 
fled so as to correlate with (N). Thus, even though the 
scatter in (N) and A is the same as in Figure |U Equa- 
tion [6^, is violated, and the clouds are inconsistent with 
M oc A by 5fi. Finally, the dynamic range in areas for 
the data in Figure[7]is artificially compressed, decreasing 
flog A- While L3 is still an appropriate fit to the data, the 
scatter about the line is larger, as suggested by Equation 

Eb. 

4. THE INTRA-CLOUD MASS- AREA 
RELATIONSHIP 

We are now in a position to address the mass-area re- 
lationship within a single cloud, and in particular the 
mass-a rea relationships presented by iKauffman n et al.1 
(I2010al). The ma in diffe rence between the analyses in 
IKauffmann et al.l (|2010aft and lLombardi etail (|2010D is 
that, in the former work, each disconnected region within 
a cloud is treated separately, while such structures are 
ag gregated in the la t ter. Th e mass-size plots presented 
in IKauffmann et al.l ()2010al ) exhibit larger scatter and 
deviations from M oc A (Figure [TJo) . Does this contra- 
dict the "universality" of cloud structure suggested by 
LAL10? Do cloud substructures exhibit substantial vari- 
ation in the column density PDF that only average out 
when entire clouds are compared to one another? 
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The formalism discussed above can be applied to sub- 
structures as well - in particular, the mass, area, and 
PDF for each substructure in Figure [TJd are still related 
via Equations [T] and [21 with the caveat that, in general, 
P{x) can vary from region to region within a cloud. How- 
ever, the lines drawn in Figure [Tp trace substructures 
within Perseus at a variety of different contour levels - 
they are an example of a Type 3 mass-area relationship 
presented in Figure [3] Consequently, they connate any 
regional change in the PDF with changes in Nq (in other 
words, the value for (N) in Equation [2] depends on both 
P(x) and N ). 

To more directly compare the internal mass-area re- 
lationship in Perseus to the inter-cloud relationship dis- 
cussed above, we show in Figure[8]the mean column den- 
sity as a function of substructure area at a variety of 
different thresholds. Each line adopts a fixed extinction 
threshold, and is thus an example of a Type 4 mass-area 
relationship. These plots possess the same features as 
those in Figure [4j in particular, changes in mean column 
density are uncorrelated with size, and criogjv <C &\ gA- 
In other words, the apparent scatter of the lines in Figure 
[Ip is due to displaying several structures defined via dif- 
ferent thresholds Nq, and the cloud substructure within 
Perseus is "universal" in the same sense as the inter-cloud 
comparison presented by LAL10. 

5. DISCUSSION AND FUTURE DIRECTIONS 

Recent studies of the mass-area relationship within 
and among molecular clouds can be understood from 
the perspective of a more fundamental quantity - 
the column density distribution function. Specifi- 
cally, the mass-area relationship serves as a diagnos- 
tic for the degree to which the PDF varies within 
and among clouds. This has a certain observa- 
tional appeal - while many different t echniques are 
used t o measure cloud ma s s and size (iSanchez et all 
2001 iRosolowskv fe Lerovl [200l iHeyer et al.l 120091: 
Kauffmann et al.ll2010at lLombardi et al.l 120101) the col- 
umn density PDF is usually measured and characterized 
in a consistent way. Thus, different observational studies 
of the PDF are more easily inter-compared than studies 
of the mass-area relationship. 

Larson's mass-area relationship is occasionally inter- 
preted as an artifact of defining cloud boundaries at fixed 
column density thresholds - the argument is that the 
mass of the region in this scenario is M ~ -/Vthrcsh x A. 
This relation is approximately valid for real clouds, given 
the fact that their column density PDFs are all similar 
and concentrated towards low column density. However, 
if column density PDFs did vary substantially from re- 
gion to region (as defined by Equation [6]) , this would 
affect the masses, and would violate L3. Larson's mass- 
area relationship is fundamentally an assertion that the 
means of column density PDFs do not correlate with the 
sizes of cloud structures. 

Volume and column density PDFs (p-PDF, A/"-PDF) 
have also received extensive theoretical attention. Log- 
normal p— PDFs are often attribute d to the action of 
isoth e rmal, supersonic turbulence ( Vaza ue z-Semadeni 
1994 ; iPadoan et al.l Il997t iPassot fc Vazquez-Semadeni 
1998). Various scenarios could map a log- normal 
p— PDFto a log-normal TV-PDF. For example, if the typ- 
ical line-of-sight cloud depth is smaller than the turbu- 



lence correlation length - and relatively uniform across 
the cloud - then the A/"-PDF would simply be a scaled 
version of the p— P DF, multiplied by the (constant) 
line-of-sight depth (I Vazquez-Semadeni fc Garcial 120011 : 
iBallesteros-Paredes et al.l 1201 1[ ). Alternatively, if cloud 
depths are large compared to the correlation length, each 
line of sight samples many independent values of the 
p— PDF. The column density along each line of sight con- 
verges to ~ (/)) x L, where L is the cloud depth and (p) 
is the mean volume density. If the distribution of cloud 
depths L follows a log-normal, so will the A/"-PDF. It 
seems plausible that real clouds lie somewhere between 
these extremes, with modulations in both cloud depth 
and volume density influencing the A/"-PDF. 

It has been suggested that the Mach number of super- 
sonic turbulence determines the the width of the p— PDF 
(jPassot fc Vazquez-S emadenl ll998T ). Since larger clouds 
exhibit larger velocity dispersions (Larson's first relation- 
ship), one might expect larger clouds to possess broader 
A/"-PDFs and, by extension, larger values of (N). Such 
a phenomenon would disobey L3, and is not seen in Fig- 
ures [4] or [H Other processes may be acting to confine 
the width of cloud c olumn density distributions, (e.g., 
magne tic turbulence. lOstriker et al.l 120011 ). iTassis et al.l 
( 2010) has also demonstrated that several different phys- 
ical processes besides supersonic turbulence can repro- 
duce the column density PDFs observed for real clouds. 

Gravity may act to create power-law excesses in the 
high column density tails of cloud PDFs. This view 
is supported observationally by the fact that non-star 
formi ng clouds tend to lack these tails (jKainulainen et all 
2009), and that the star formation rate correlates with 
the integrated ma ss of the hig h-column density regions 
(|Lada et al.l l2010l ). However. lAlves. Lombardi fc Lad"al 
(|2012[ ) suggest that these tails may also be consistent 
with the superposi t ion of several cloud components, and 
iKainulainen et all (|2011l ) have subsequently suggested 
that the high-column density tails are confined by pres- 
sure instead of gravity. Substructures within a cloud 
may well be dominated by different physical processes, 
inducing regional changes in the PDF. Figure suggests 
that the mean of the PDF does not change dramatically 
within Perseus. 

While the mass-area relationship can be viewed as 
a probe of the column density PDF and its variation, 
it contains only limited information about the PDF 
(namely, its mean). Given the amount of theoretical con- 
sideration of the A/"-PDF, more carefully characterizing 
this distribution may provide more insight into what in- 
fluences cloud structure. A rigorous observational study 
of the A/"-PDF and its variations, however, must account 
for two factors: 

1. Extinction measurements are still fairly noisy, and 
observed column density distributions convol ve the 
true P DF with the measurement error kernel (jKellvl 

EqUJ). 

2. Clouds exhibit spatial correlations, and so neigh- 
boring measurements of the column density do not 
constitute independent samples of the underlying 
A/"-PDF. Fitting procedures which do not account 
for this will over-estimate the effective number of 
samples, and under-estimate uncertainties in the fit 
to the PDF. 
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TABLE 1 

Literature Fits to (M/Mq) = 7(A/pc 2 )' 3 
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P 


Reference 


Tracer 


Type a 


240 


0.95 


Larson (1981) 


13 CO, 12 CO, NH 3 


1 


170 


1 


Solomon et al. (1987) 


12 CO 


1 


30 


1-1.2 


Falgarone et al. (1992) 


12 CO, 13 CO 


1 


50-500 


1.2-1.8 


Elmcgreen & Falgarone (1996) 


12 CO, 13 co, c 18 o 


1 




0.91 


Sanchez et al. (20051 


13 co 


3 


150 


1.3 


Lada et al. (2008) 


Extinction 


1 


42 


1 


Hever et al. (2009) 


13 CO 


1 


228 


1.36 


Roman-Duval et al (2010) 


12 CO, 13 CO 


1 


70 


0.5—1 


Kauffmann et al. (_2010a) 


Extinction, mm-continuum 


3 


41-380 


0.99 - 1.01 


Lombardi et al. (2010) 


Extinction 


2 


150 


1.3 


Roman-Zuniga et al. (2010) 


Extinction 


1 



a Mass-area comparison type, as defined by Figure [3] 



To date, most characterizations of cloud PDFs have ig- 
nored these issues. However, they must be addressed to 
properly constrain the shape and spatial variation of the 
PDF. 

6. CONCLUSION 

We have directly compared the mass-area relationship 
within and among clouds, interpreting this relationship 
from the perspective of the column density distribution 
and its spatial variation. For structures defined via a 
fixed extinction contour, the mean of the column density 
PDF in each region varies less than, and is uncorrelated 
with, the region-to-region dispersion in area. This nat- 
urally yields the relationship M oc A, and thus suggests 



that Larson's third scaling relationship holds both within 
and among clouds. 

The structure of nearby molecular clouds is "univer- 
sal" in the sense that the dispersion of the PDF mean 
is negligible compared to the dispersion of cloud areas. 
However, cloud structure may still vary in interesting 
ways that do not greatly affect the mean of the A^-PDF. 
We suggest that studying the shape and variation of the 
A/"-PDF directly may provide more insight into the uni- 
versality of cloud structure. 

This material is based upon work supported by the 
National Science Foundation under Grant No. AST- 
0908159. CRZ acknowledges support from program 
UNAM-DGAPA-PAPIIT IA101812-2. 
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Left: The inter-cloud mass-radius relationship, as published by Lornbardi et al. (2010). The two lines correspond to cloud 
defined at to different extinction contours Aq(= Ak ), The data ob e y M oc B? to excellent approximation. Right: The intra- 
-radius relationship for Perseus, as presented by Kauffmann et al. (2010a). Each line corresponds to a single substructure, 

different extinction contours. Different substructures appear to have significantly different relationships. Both figures define 




Fig. 2. — Schematic comparison of how cloud masses/areas are m easured by [L ornbardi et al. (2010) and Kauffmann et al. (2010a). Each 
closed c ontour has a well -defined mass and area, as described above. Kauffmann et al. (2010a) treat each closed contour separately (circles), 
whereas Lornbardi c i a j (2010) add up the mass and area of every region defined by a given intensity threshold (squares). 
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Fig. 3. — The different possible mass-area relationships to measure. As in Figure[2] each point corresponds to a mass and size measurement 
for a single closed contour in a cloud. The color of the point represents the contour value, and the shape represents the cloud. In this paper 
we focus on the relationships 2 and 4; that is, our comparisons of structures within and betwe en clouds always adopt a fixed extinction 
threshold. The numbering scheme has been chosen to match the linewidth-size types defined by Goodman ct al. (1998). 




Fig. 4. — The connection between the column density PDF and mass-area relationship for the data presented inXombardi et al. (2010). 
Left: the column density PDF for each cloud. Right: the mean column density of pixels with extinctions higher than the threshold defined 
by the solid line on the left plot, as a function of cloud area. Bottom: The mass-area relationship. 
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Fig. 6. — The same as Figure [4] but with the cloud areas shuffled so as to correlate with mean column density. These clouds are 
inconsistent with L3 at the 5cr level, even though the spread in mean column density and area is the same as in Figure [4] 
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Fig. 7. — The same as Figure [4] but with the dynamic range in areas artificially compressed. L3 is still an appropriate model, but the 
statistical errors are larger. 
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Fig. 8. — The mean column density and area for substructures in the Perseus cloud, at three different column density thresholds. As in 
the inter-cloud relation, the spread in mean column density is small compared to the spread in area, and uncorrclated with area. 



